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Abstract 

Semimicroscopic optical potentials for low energy proton reactions in 
mass 90 — 100 region have been obtained by folding the density dependent 
M3Y interaction with relativistic mean field densities. Certain parame- 
ters in the potential have been deduced by comparing calculated results 
with the data for elastic scattering. Low energy proton reactions in this 
mass region have been studied in the formalism with success. Rates of 
important astrophysical reaction in the mass region have been calculated. 

The p-process is a common term given to the astrophysical reactions, which 
are involved in synthesis of heavy elements but do not correspond to the r- or 
the s-processes. It includes reactions such as proton capture, charge exchange 
and photodisintegration. The p-process is known to be important for production 
of certain so called p-nuclei, which are beyond the ambit of the slow and fast 
neutron reactions. A p-network involves typically two thousand nuclei, and in- 
corporates, again typically, twenty thousand reactions and decays. More details 
may be found in standard text books [for example Illiadis [1] ] and reviews P]- 

One obvious problem in studying the p-process is that many of the involved 
nuclei have very short life times and are not available in our terrestrial laborato- 
ries for experiment. Though radioactive ion beams have opened a new vista, we 
are still far away from having the reaction rates at astrophysical energies for all 
the main reactions involved in the p-process. Thus, theoretical calculations for 
such rates remain very important for the p-process. For example, Rapp et al. 
have identified a number of reactions, which are very important in p-process[3 . 
In mass 90 -100 region, the list includes the photodisintegration reactions emit- 
ting protons and leading to the products 91 Nb, 95 Tc, and "Rh. The rates for 
these reactions could be obtained from their inverse, i.e. (p,j) reactions, had 
the above nuclei been easily available in laboratories. 

There have been numerous theoretical calculations of astrophysical rates [1] 
employing various models. However, very often in literature, these theoretical 
rates are varied by factors ranging from ten to hundred to study their effects [S]- 
In some earlier works, we calculated the cross sections of various low energy 
proton reactions some of which are involved in the rapid proton processes in 
mass 60-80 regionjBl |S]- The semimicroscopic optical model was employed 
for calculation of cross sections using densities from theoretical mean field cal- 
culations. It is our aim to fix the various parameters and prescriptions in our 
procedure by fitting available low energy cross sections for various reactions in 
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a mass region and calculate the rates for various reactions involving protons, 
which arc important in nucleosynthesis. Thus, a more stringent restriction may 
be imposed on the variation of rates. Some consequences of the above approach 
in rapid proton process have already been analyzed [91 110). 

The code TALYS1.4[TT] has been used to calculate cross sections and rates 
in the Hauser-Feshbach formalism. In our earlier works, it was concluded that 
the Hartree-Fock-Bogoliubov level densities, calculated in TALYS by Hilaire[T2"] 
and the El gamma strength functions, calculated in the same approach, fit the 
results in mass 60-80 region. In the present calculation we employ these values 
to extend our calculations to mass 90 region. 

The method followed in the present procedure has been detailed in our earlier 
publications [3 [7] and is not discussed here. The FSU Gold[T3] Lagrangian 
density is employed to calculate the nuclear density. Since we need the density 
as a function of radius, the calculation is performed in the co-ordinate space. 
We employ spherical approximation, as most the nuclei under study are near 
closed shells for both protons and neutrons and are not strongly deformed. The 
effective interaction DDM3Y[T3], derived from nuclear matter calculation, has 
been folded with the nuclear densities to obtain the semi-microscopic optical 
model potentials in the Local Density Approximation. 

Since the nuclear density is an important factor in the present formalism, we 
have studied the charge radii values. The charge radius is the first order moment 
of the charge distribution. In Table[TJ we compare our results for the charge radii 
( r ch) with measurements for those nuclei in this mass region, which have been 
involved in the reactions studied later in this work and for which experimental 
radius information are available. Charge densities have been obtained by folding 
point proton densities with a Gaussian form factor to incorporate the effect of 
the finite size of the proton as in our previous work [8]. It is clear that the charge 
radii are reasonably well produced in our calculation. 

We could not find direct experimental values for charge densities. Hence, 
we have employed the Fourier-Bessel coefficients for densities extracted from 
electron scattering experiments in de Vries et al. [16] to get the charge densities 
and plotted two examples in Fig. [T] One can see that the theoretical results 
reasonably agree with experiments. However, the absence of any information 
on error prevents us from reaching a firm conclusion. 

As a first test of the optical model potential, we have looked at elastic 
proton scattering at low energies. Elastic scattering involves the same incoming 
and outgoing channel for the optical model and may be taken to provide the 
simplest test to constrain various parameters involved in the calculation. The 
proton energy relevant to a typical p-process temperature of 1—3 GK for nuclei 
in this mass region lies between 1—4 MeV. However, scattering experiments are 
very difficult at such low energies, as the cross sections are extremely small, and 
hence no experimental data are available. We have compared the cross sections 
at the lowest energies available in literature with theoretical results. 

In Figs. [HandGIl we present the results of some of our calculations in Zr and 
Mo isotopes, respectively, along with the corresponding experimental results. 
Experimental values are respectively from Refs. [TT1 [TS1 US] for 90 > 91 > 92 Zr, and 
from Ref. [3D] for Mo isotopes. To fit the experimental data, the folded DDM3Y 
potential has been multiplied by factors of 0.81 and 0.15 to obtain the real and 
imaginary parts of the optical potential, respectively. Throughout the rest of 
the work, we use these two factors to obtain the potential. We emphasize that 
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better fits for individual reactions are possible by varying different parameters. 
But if the present calculation has to be extended to unknown mass region, this 
approach is clearly inadequate. Therefore, we have refrained from fitting indi- 
vidual reactions. In our previous work[71IH], we used a different normalization 
which is in good agreement with experimental values in a wide mass region(A 
« 60-88). But beyond that region, same set of parameters are unable to fit 
the experimental data for p- nuclei [3] and therefore, we choose the above set of 
parameters. Though, there are no sharp boundaries for a mass region, but for 
simplicity, we choose it in such a way that a single set of parameters can fit the 
entire mass region. In present work, we have chosen the mass region A«89-100. 

It is clear from Figs. [2] and [3] that the DDM3Y interaction can describe the 
data well. In fact, we have found that as one goes to lower energies, the quality of 
agreement tends to improve. Thus at energies relevant to astrophysical interest, 
we can expect the present method to provide a good description. 

The same formalism has been used to study the low energy (p, 7) reactions in 
a number of nuclei in this mass region. As the cross section varies very rapidly 
at low energies, it is more convenient to present the S-factor values. In Figs. 
U[71 calculated values are compared with experimental results. Next, we very 
briefly discuss our results. 

For 89 Y, the experimental values are from Tsagari et al. [H]. For 96 Zr, the 
results are from Chloupek et al. [22] , though there seems to be certain error in 
the values in that reference. The numerical values presented there are larger by 
a factor of 10 3 that the values presented in Fig. 9 of the reference. The latter 
values appear to be correct to the present authors and and are indicated in Fig. 
IH The data for Mo and Ru isotopes are from Ref [533 < respectively. 

For 96 Zr target, there are very few experimental points within the energy 
range important for astrophysical reactions. In 89 Y, the trend of the experi- 
mental values has been correctly reproduced. None of the experimental values 
differs by a more than a factor of two. The last two comments are generally valid 
for almost all the other reactions. One important exception is the 98 Ru(p, 7) re- 
action, where the measured cross section systematically increases with decrease 
in energy below 2 MeV proton energy compared to the calculated values and 
becomes larger by more than one order of magnitude around 1.6 MeV. It has 
not been possible to explain such a large increase, which is absent in all low 
energy (p, 7) reactions in this mass region for which data are available. In fact, 
the calculation carried out in Ref. [24] , where the experimental values have been 
published, is also unable to explain such a sudden increase. 

We next look for other low energy reactions involving proton projectile. The 
only reaction for which we have been able to find substantial amount of data in 
the domain of astrophysical energies is the 93 Nb(p, n) reaction[25j. Our results 
have been presented in Fig. [51 One can see that our calculation gives an 
excellent description of the experimental trends. However, one should also note 
that the data are rather old and have either very large errors or no quoted error 
value. 

For the sake of completeness, in Fig. Owe compare the rates of (7,p) reaction 
from present calculation with rates from NON-SMOKER^] calculation for 92 Mo 
and 96 Ru. One can see that the present calculation is very similar to the NON- 
SMOKER values. Therefore, it is expected that all the results can also be 
reproduced with commonly used NON-SMOKER rates. 

From the above discussion, it is possible to conclude that the low energy 
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reaction cross sections are, except in one case, reasonably reproduced in the 
above approach. The success in this calculation has enabled us to calculate the 
astrophysical rates for the reactions, identified as important by Rapp et at [3], 
in mass 90-100 region. They are presented in Table [2] 

To summarize, relativistic mean field calculation has been performed in nu- 
clei between mass 90 and 100 to obtain the density profiles. They, in turn, 
have been folded with the density dependent M3Y interaction to obtain the 
semimicroscopic optical potential. Parameters in the potential have been fixed 
by comparing with low energy proton scattering. Available experimental infor- 
mation on low energy proton reactions has been compared with theory. Rates of 
important astrophysical reaction in the mass region have also been calculated. 
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Table 1: Experimental charge radii values compared with calculated results for 
the nuclei involved in low energy proton reactions. The experimental values are 
from the compilation by Angeli[15]. 





Teh 


(fin) 




r c h 


(fm) 




Thco. 


Exp. 




Theo. 


Exp. 




4.274 


4.242 


y4 Mo 


4.366 


4.352 


90 Zr 


4.297 


4.270 


95 Mo 


4.376 


4.362 


92 Zr 


4.317 


4.306 


96 Mo 


4.387 


4.384 


94 Zr 


4.335 


4.331 


98 Mo 


4.407 


4.409 


96 Zr 


4.356 


4.350 


96 Ru 


4.409 


4.393 


92 Mo 


4.344 


4.316 


98 Ru 


4.431 


4.409 



Table 2: Rates in cm 3 mole 1 sec 1 for selected (7,p) reactions of astrophysical 



importance. 


T (GK) 






Target 






92 Mo 




96 Ru 


ioo pd 


1.5 


3.45x10" 


-15 


3.67xl0~ ib 


2.69X10" 14 


2.0 


3.34x10" 


-07 


3.64xl0-° 7 


1.52xl0-° 6 


2.5 


2.89x10" 


-02 


3.31xl0~ 02 


9.85xl0~ 02 


3.0 


6.77x10" 


-01 


7.95 xl0+ 01 


1.85xl0+ 02 


3.5 


1.92x10" 


-04 


2.23xl0+ 04 


4.30xl0+ 04 


4.0 


1.40x10" 


-06 


1.57xl0+ 06 


2.60xl0+ 06 


5.0 


5.87xl0+ 08 


5.70xl0+ 08 


7.86xl0+ 08 
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Figure 1: Comparison of charge density obtained from Fouricr-Bcssel analysis of 
experimental electron scattering data (solid line) and calculated in the present 
work (dashed line) for (a) 90 Zr and (b) 94 Mo respectively. 
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Figure 2: Experimental and calculated cross sections for elastic proton scattering 
in Zr isotopes. For A = 90, 91 and 92, the proton energies are 9.7 MeV, 14.8 
MeV, and 14.25 MeV, respectively. The cross sections for 91 > 92 Zr have been 
multiplied by a factor of 100 and 1000, respectively. 
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Figure 3: Experimental and calculated cross sections for elastic proton scattering 
in Mo isotopes at 15 MeV proton energy. The cross sections for A = 92, 94, 
and 96 have been multiplied by factors of 10, 100, and 1000, respectively. 
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Figure 4: Experimental and calculated S-factors for (p, 7) reactions in 89 Y 
and 96 Zr, respectively. The solid (dashed) line indicates calculated results for 
89 Y( 96 Zr). 
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Figure 5: Experimental and calculated S-factors for 92 ' 94 Mo(p, 7) reactions. 
Results for 94 Mo have been multiplied by 10. The solid (dashed) line indicates 
calculated results for 92 Mo( 94 Mo). 
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Figure 6: Experimental and calculated S-factors for 95 ' 98 Mo(p,7) reactions. 
Results for 98 Mo have been multiplied by 10. The solid (dashed) line indicates 
calculated results for 95 Mo( 98 Mo). 
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Figure 7: Experimental and calculated S-factors for 96 ' 98 ' 99 Ru(p, 7) reactions. 
Results for 98 ' 99 Ru have been multiplied by 10 and 100, respectively. The solid, 
dashed and dotted lines indicate results for 96 > 98 >"R U) respectively. 




Figure 8: Experimental and calculated S-factors for the 93 Nb(p, n) reaction. 
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Figure 9: Comparison of rates (cm 3 mol^ 1 sec -1 ) for (j,p) reactions from 
present calculation (solid line) and NON-SMOKER[4 calculation (dotted line) 
for (a) 92 Mo and (b) 96 Ru respectively. 



12 



